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We investigate the clustering properties of a dynamical dark energy component. In a cosmic mix of 
a pressureless fluid and a light scalar field, we follow the linear evolution of spherical matter pertur- 
bations. We find that the scalar field tends to form underdensities in response to the gravitationally 
collapsing matter. We thoroughly investigate these voids for a variety of initial conditions, explain 
the physics behind their formation and consider possible observational implications. Detection of 
dark energy voids will clearly rule out the cosmological constant as the main source of the present 
acceleration. 
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I. INTRODUCTION 

Direct [l[ and indirect Q evidence of the present ac- 
celeration of the universe is accumulating. Nonetheless, 
the source of the accelerated expansion is as elusive as 
ever. In the context of general relativity, one needs to 
hypothesize a dark energy source with negative pressure 
to accommodate the acceleration. Alternative ideas in- 
clude various modifications to gravity. Little is known 
about the dark energy except that its pressure is nega- 
tive, and that it accounts for about 70% of the critical 
density of the universe. Whether the dark energy is the 
cosmological constant or a dynamically evolving source 
of energy is a fundamental unanswered question. While 
the cosmological constant has a fixed ratio of pressure to 
energy density, w = p/p = —1, dynamical dark energy 
(DDE) will in general have a varying equation of state 
(EOS), w{z). Observing a deviation from —1 or a time- 
evolution in the EOS will be decisive evidence in favor of 
the existence of DDE. However, there are known degen- 
eracies [|[ which make this task extremely difficult, un- 
less the deviation from a cosmological constant is strong. 
The current observational limit on the EOS of dark en- 
ergy is roughly w « — 1 ± 0.1 at the la level [lj, which is 
consistent with a cosmological constant. Future experi- 
ments hold out the possibility of pinning down this limit 
by maybe a factor of 10. For a recent review, see 

The effect of the cosmological constant and DDE on 
the expansion rate can be identical, and there is a need 
for probes that go beyond the background to distinguish 
between the two. An example of such a probe is structure 
formation. There are numerous works exploring the for- 
mation of structure in the presence of homogeneous dark 
energy [5j . An exciting and somewhat controversial pos- 
sible difference between the cosmological constant and 
DDE is their clustering behavior. While the cosmologi- 
cal constant is exactly homogeneous on all scales, DDE 
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is expected to be not perfectly homogeneous [6|, and the 
implications of this on the CMB are well known ,7j] . How- 
ever, it is usually assumed that the clustering of DDE is 
negligible on scales less than 100 Mpc. Whether small 
perturbations in DDE can be neglected is debatable, and 
a deeper understanding of the DDE inhomogeneous dy- 
namics is clearly needed. 

Several recent works have explored the consequences 
of DDE clustering on scales shorter than 100 Mpc. Some 
have adopted a phenomenological ajpproach, parameteriz- 
ing the clustering degree of DDE [8| . These works point 
out potential observables of DDE clustering, justifying 
further investigation. Works which attempt a more fun- 
damental treatment are mostly in the context of coupled 
dark energy [§] , or other non-trivial models of DDE [lol ] , 
as clustering is most probable in such theories. However, 
less attention has been given to the clustering in simpler 
models of DDE. 

In the present work, we aim to further explore the inho- 
mogeneous behavior of DDE. Our approach is straight- 
forward: starting with a gravitational action which in- 
cludes matter and DDE, we numerically follow the lin- 
ear evolution of spherical perturbations of matter and 
the DDE response to these perturbations. For the sake 
of simplicity, our model for the DDE is a light scalar 
field, which is not explicitly coupled to the matter den- 
sity. As the only coupling between the DDE and the 
matter is gravitational, our results are conservative in 
the sense that any model more complicated will exhibit 
stronger DDE perturbations than shown here. The strik- 
ing feature that emerges from our calculation is that in 
the vicinity of collapsing matter, the DDE develops a 
spatial profile and tends to form voids. The mechanism 
that allows the void to form is that although initially the 
field's evolution is friction dominated due to the cosmic 
expansion, the collapse of matter slows down the local 
expansion. This allows the field to locally roll down and 
lose energy, creating the void. The presence of the matter 
perturbation is necessary to trigger this mechanism. 

The plan of the paper is as follows: in section [H] we 
describe our model in detail. In section llHl we present our 
results. We discuss our results and state our conclusions 
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in sections IIVI and [V] respectively. 



ordinate, except for V '(</>) 
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II. THE MODEL 

We are interested in spherical perturbations around a 
flat FRW universe. The most general line element in 
comoving coordinates is then 



ds 2 



dt 2 -U(t,r)dr 2 -V(t,r)(d6 2 + sin 2 9dip 2 ) , (1) 



where U{t,r) and V(t,r) are general functions [Tl|. 

We take a cosmic mix of non-relativistic matter and a 
DDE component as the energy source. The matter com- 
ponent is described by a perfect and pressureless fluid, 
with an energy-momentum tensor given by 



T^(m) = diag(>, 0,0,0) 



(2) 



where p is the energy density of matter. 

We model the DDE with a classical scalar field cj) with 
a Lagrangian C given by 



C = \ {d^f - V{4>) 



and an energy-momentum tensor given by 
The EOS of the DDE w is defined as 
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These are subject to the following constraint equations: 
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Linearization 



with the energy density p^ and the pressure are 
read off the energy momentum tensor, Tbo(^>) and 
— g l ^Tij(4>)/3 respectively. 

It is convenient to rewrite Einstein's equations in the 
following way, 



(6) 



where i? M „ is the Reimann tensor, and K = 8wG. As 
there is no explicit interaction between the matter and 
the DDE, energy conservation applies to each separately, 



V M T^(m) = 0, V M r^(0) = O 



(7) 



The time-evolution of the system is given by the fol- 
lowing equations (where dots denote time-derivatives and 
primes denote derivatives with respect to the radial co- 



We now proceed to separate our variables to a homoge- 
neous background and a time and space-dependent per- 
turbation, which we will then linearize. Working in the 
synchronous gauge fl2T | , we redefine the metric functions 
U and V as follows: 



U(t,r) = a(t) 2 e 2 ^ 
V(t,r) = r 2 a(t) 2 e 2 ^ 



(14) 



Here a(t) is the scale factor of the spatially homogenous 
and flat background, and C(t, r) and i/j(t, r) are the devi- 
ations. We introduce a perturbation around a homoge- 
neous background also in the matter and the DDE, 

p(t,r) = p(t)+5p(t,r) 
cj>(t,r) = <j)(t) + S^{t,r) 
V{(t> + 84>) = V(cj))+SV((t>,Scj)) . 
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The zeroth order of equations j8]) 
3H 2 -K p + V+-^ 



13]) gives 

= (15) 



H + 3H 2 - K + V^j = (16) 

p + iHp = (17) 
4> + 3H4> + V = , (18) 

where H = a/a is the Hubble function. 

To linear order, the evolution equations (f8j)- ([TTj) give 

•• • -2 (C! 2ib' 
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and the constraint equations (fl2l) - p3]l reduce to 

2H( + + + ^ " tp 

a \r r r r 

-K Up + 5V + <M0) = 
- (i> + C) - + 2ib' = 

Combining equations ([!§]) , ([2"0]) and ([2"3"]) gives 

C + 2$) + 2H U + 2ijj 

+K(Sp- l )\' + 2,,,),,) = 



(23) 
(24) 



(25) 



The only combination which is relevant to the equa- 
tions of motion (j2"Tj) and (|2"2"|) is x = C + 2j/>. Compar- 
ing ([17]) and ([2"T]) it is clear that \ can be thought of as 
3SH, and therefore characterizes the spatial profile of the 
Hubble function. At the cost of losing some information 
about the metric, we can reduce the number of our vari- 
ables and equations from 4 to 3 by solving for x instead 
of for C and ip. Equations pi]). ([22]) and ([25]) yield 



(26) 
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+ 3H6j> + SV + x4>- 4V 2 (5ct>) = (27) 



X + 2ff X + K ( -<Jp - SV 
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Finally, by Fourier-transforming 8p(t,r), S(f>(t,r) and 
x(t, r) into 5pk(t,k), S<bk{t,k) and Xfc(^^) respectively, 
equations ([26]) -([28 ]) can be written as a set of ordinary 
differential equations: 

6p k + 3HS Pk + p X k = (29) 



6,. + 3H5<p k +[V" + 



k + 4>Xk=Q (30) 



Xk + 2H Xk + K i^-Spk - V'5<Pk + 20 k 5<Pk) = , (31) 

where we have used the fact that to linear order, SV — 
V'5cj> and SV = V"5<j>. 



B. Potential 

Observationally distinguishing between various poten- 
tials of DDE is a formidable task [13], and a careful 
analysis of the growth of structure in various potentials 
might prove a useful tool. Our present goal though is to 
trace generic properties of DDE. Accordingly, we choose 
to work with a simple mass potential, 



V{4>) = -m 2 <\> 2 



(32) 



We take the mass scale comparable to the present Hubble 
scale, m<p/Ho ~ 1. The light mass assures a slow roll be- 
havior, which will provide accelerated cosmic expansion. 
Unless noted otherwise, the figures presented here refer 
to this mass potential. 

In order to verify the generality of our results, we re- 
peated the analysis for a more realistic potential - the 
double exponential 



v{4>) = v (t 



(33) 



with a = 20.1 and f3 — 0.5. As we later show, the result- 
ing behavior for the two potentials is qualitatively the 



C. Initial conditions 

We want to study how the DDE reacts to the clustering 
of matter. Thus our initial conditions are of perturbed 
matter and homogeneous DDE. The matter perturbation 
is taken as a spherical gaussian, 

S m (ti,r) ee Sp m (ti,r)/p m (ti) = Aexp (~r 2 /a 2 ) , (34) 

and the DDE is taken to be initially homogeneous, 8<fi — 
0. A non-homogeneous evolution for the dark energy 
is nonetheless allowed. To ensure that the matter per- 
turbation has no peculiar velocity, we choose the initial 
condition for the metric variable as x — 0. From equa- 
tion (1261). we find that this amounts to the statement 
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that initially, 8p oc a 



the matter particles mak- 



ing up the perturbation are being simply carried along 
with the Hubble expansion (or in other words, they are 
at rest in a comoving frame). However, we have veri- 
fied that choosing different initial values of \ does not 
qualitatively affect our results. 

The numerical calculation begins at a rcdshift of z = 35 
with a perturbation amplitude of A = 0.1, and is run to 
the present, z = 0. We focus on relatively short-scale 
perturbations, aHi ~ 0.01. With this choice of redshift 
span, we aim to explore the full linear range of the matter 
perturbation, while also allowing for a period of DDE 
dominance. 

The initial values of the background variables pi, 4>i 
and Hi are chosen such that their present values (denoted 
with a subscript) converge to fl m ,o = 0.3, O^o = 0.7, 
and a normalized Hubble value, Hq = 1. This is done 
with the help of a root-finding algorithm. We take <j>i = 0, 
which means that the initial state of the scalar field is 
homogeneous and with equation of state w = — 1, similar 
to the cosmological constant. 

We now have the layout to numerically evolve equa- 
tions (p^ - lpTf]) . The solutions for Spk, 5<t>k and Xfe are 
then Fourier transformed back to real space. 



III. RESULTS 
A. Density contrasts 

The results shown in this section are for the mass po- 
tential, equation (|32[) . unless specifically noted otherwise. 

The numerical run begins at a redshift of z = 35, when 
the DDE is subdominant. The DDE remains subdomi- 
nant through most of the growth time of the perturbation 
as well, and accordingly, we expect the matter density 
contrast 8 m = Sp m /p m to grow as the scale factor a. 
This is indeed confirmed by our results. 

The most striking result that emerges from our cal- 
culation is that in the vicinity of collapsing matter, the 
DDE tends to form voids, or growing regions of under- 
density. The anti-correlation between the perturbations 
of the matter and the DDE was noted in Q and in [l5[ 
for particular cases. Figure [T] plots the DDE density con- 
trast ($0 = Sp^/p^ at the center of the matter perturba- 
tion, r = 0, against the redshift (1 + z). We find that 
the amplitude of the perturbation grows sharply at late 
times. 

Figure [2] shows the growth of the absolute value of the 
DDE perturbation with redshift on a logarithmic scale. 
This plot reveals another interesting effect: the initial re- 
sponse of the DDE to the gravitationally collapsing mat- 
ter is a very weak tendency to collapse. The collapsing 
phase, however, is extremely short-lived (an O[10 _1 ] frac- 
tion of the total time of the run) . The kink in the figure 
is the cross-over from a positive to negative perturbation. 

We next look at the spatial profiles of the perturba- 
tions. Figure [3] shows the profiles of both the S m and 8$ 
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FIG. 1: The evolution of the DDE overdensity 8$ at the cen- 
ter of the matter perturbation, r = 0, with redshift (1 + z). 
The scale of the perturbation is aHi — 0.01, and the mass 
is m/Ho = 1. Initially homogenous, the DDE develops an 
underdensity at late times in response to the matter pertur- 
bation. 
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FIG. 2: Same as Figure [T] with the y-axis on a logarithmic 
scale. The DDE tends to cluster initially, but eventually forms 
a void. 
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three different redshifts. Solid lines denote the DDE profiles, 
and dotted lines denote the matter profiles. 



3: Logarithmic profiles of the matter density contrast 
10 1 5m | and the dark energy density contrast log 10 l^l, at 
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at an early stage of the run (z = 27.5), at an intermedi- 
ate redshift (z = 0.8) and at the final stage [z = 0). The 
x-axis shows the physical scale as a fraction of the hori- 
zon size. The amplitude of the perturbations is shown 
on a logarithmic scale, and it is worth noting the change 
in gap between the two scales. The growth rate of 5$ is 
significantly faster than that of 5 m , so that the two ampli- 
tudes are almost comparable at late times. This suggests 
that a calculation of the non-linear regime might reveal 
interesting behavior of the DDE. 

Figures H] and [5] show how the growth of the DDE per- 
turbation is affected by the initial width of the matter 
perturbation. The figure confirms the sensitivity of DDE 
perturbations to the scale, and as expected, shorter scales 
exhibit a suppressed behavior. Nonetheless, a possibly 
significant amplitude can be found on relevant scales. 
For instance, our local supercluster has diameter of or- 
der 0.01-ffg" 1 and an overdensity of about S m w 1.2. Our 
runs of aHi = 0.01 end with an overdensity of S m w 1.6, 
but roughly they can be taken as a measure of what we 
should expect on these scales. Also, one must remember 
that this analysis only captures the physics of the linear 
regime. Given the sharp increase in the void amplitude 
at later times, it is not unreasonable to expect interesting 
effects in the strongly non-linear regime. 

Figures [5] and [7] examine the sensitivity to the field's 
mass. As expected, increasing the mass of the scalar field 
causes 5$ to grow stronger and at increasingly earlier 
redshifts. One can also see that if the mass is heavy 
enough (an order of magnitude larger than the Hubble 
mass), the field will have had enough time to enter a 
period of rapid oscillations, and will effectively behave as 
regular matter. 



B. Equation of state 

We would now like to focus on the local behavior of 
the equation of state (EOS) of the scalar field. Let us 
define wq as the background homogeneous EOS, Sw as 
the leading order correction to wq and let the first-order 
corrected w be w\ = wq + Sw: 
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(36) 
(37) 

-SV, 



where to first order Sp = 4>5<f) + 5V and Sp = 
and we have suppressed the (f> subscripts. 

Figure [5] shows how w% increases with time and 
becomes non-homogenous. The EOS at the perturbation 
is less negative than the background, but for this choice 
of mass it is still negative enough to behave as dark 
energy. 




FIG. 4: DDE density contrast 6$ at the center of the matter 
perturbation r = as a function of the redshift (1 + z) for 
fixed mass m/Ho = 1 and different initial matter perturba- 
tions' widths. The larger the initial matter perturbation, the 
stronger is the void. The curves of aHi = 0.01 (dashed) and 
0.1 (solid) almost overlap. The figure zooms on late times, 
z < 3. 
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FIG. 5: Same as Figure [4] with the y-axis on a logarithmic 
scale. The shorter scales start evolving at later times than 
the longer scales, but their evolution is faster. The curves of 
aHi = 0.01 (dashed) and 0.1 (solid) almost overlap. 



To quantify the extent of the inhomogeneity in w we 
define A w to characterize the % deviation of the local w 
from the background value, 
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(38) 



The evolution of the A w spatial profile is shown in Figure 

El 



IV. DISCUSSION 



Void formation 



We now aim to intuitively explain why the DDE per- 
turbation initially increases and then sharply drops in 
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FIG. 6: The DDE density contrast 5^ at the center of the 
matter perturbation, r = 0, against redshift (1 + z) for a -Hi = 
.01 and three different masses. The figure zooms on late times, 
z < 7. 
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FIG. 7: Same as Figure [6] with the y-axis on a logarithmic 
scale. The perturbation's is extremely sensitive to the mass 
scale. 

the presence of the matter perturbation. Initially, both 
the scalar field and the Hubble function are homogenous. 
The matter perturbation introduces an inhomogeneity in 
the Hubble function, as regions with matter overdensity 
expand slower due to the increased gravitational pull. 
As a result, H acquires a spatial profile which evolves in 
time along with the matter overdensity, due to the p\ 
term in equation (f!?rJ)) . This is illustrated in Figure [TOl 
which plots the spatial profile of the % deviation of H 
from the homogeneous background value, at three differ- 
ent redshifts. Regions which have a lower local value of 
H offer less Hubble damping to the scalar field, as the ct>x 
contribution in equation (|27[) is negative. Therefore, in 
these regions the scalar field accelerates down its poten- 
tial slightly faster, the local 4> + 8<f> has a bigger absolute 
value than the acceleration of the background, (j). Thus 
the matter perturbation imparts a local "downhill kick" 
to the field, and the strength of the kick depends on the 
magnitude of the matter perturbation at that point. The 
presence of matter is essential to trigger this mechanism. 
A similar conclusion was reached in [lfj. 
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FIG. 8: Plot of wi vs r for three different redshifts. 
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FIG. 9: Plot of % deviation in w vs r at three different red- 
shifts. 



Since the matter perturbation is gaussian, we can ex- 
pect the profile of 8<j> to be gaussian as well, which is 
confirmed in figure 1 1 1 1 The acceleration of the field per- 
turbation 8(f) is zero initially (z = 35), but it quickly 
takes on a gaussian profile, which initially grows, and 
then shrinks at later times. The important point to note 
in figure QT] is that the acceleration weakens, but doesn't 
change its sign throughout the evolution. 

The acceleration profile leads to spatial variations in 
the energy density of the DDE. For the mass potential, 
the linear order of the density perturbation is 8ptf, = 
<t>8<j) + m 2 if)8(l), where the first term represents the local 
variation in kinetic energy (KE) and the second term 
represents the local variation in potential energy (PE). 
Whether the initial value of the scalar field was shifted 
to the right or left of its minimum, <f> and 8<f> will be 
of the same sign, assuring a positive correction to the 
KE. On the other hand, <f> and 8<f> will be of opposite 
signs, assuring the PE correction will be negative. The 
total correction to the energy density then will be pos- 
itive when the KE dominates, creating an overdensity, 
and negative when the PE correction dominates, creat- 
ing a void. The fact that 8<f> weakens but doesn't change 
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FIG. 10: Percentage variation of the local Hubble parameter 
at three different redshifts 



FIG. 12: The equivalent of figure [3] with the double exponen- 
tial potential, equation (|33[) . 
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FIG. 11: Profile of S(j> at three different redshifts 

its sign assures that 5<f> approaches a constant value, but 
6<fi keeps growing. The PE contribution m 2 c/><5</> soon be- 
comes dominant over the KE correction, 4>8<f). As soon 
this happens, a void is created. 

Similar reasoning should apply to other slow-roll poten- 
tials. 



FIG. 13: The equivalent of figure|8] with the double exponen- 
tial potential, equation (|33[) . 

growth of its amplitude is significantly faster than that 
of the matter perturbation, motivating a further investi- 
gation into the non linear regime. 



V. CONCLUSIONS 



B. Generality 

The results we have presented so far were for the mass 
potential, equation (|32p . This is not a very attractive 
model from the theoretical point of view, as the choice of 
the mass scale and initial conditions are fine tuned. 

While the general problem of fine tuning has not yet 
been resolved, the issue of the initial conditions is allevi- 
ated in tracking potentials such as equation (|33[) . To ver- 
ify that our results are not unique to the mass potential, 
we have repeated the analysis for the double exponen- 
tial potential. We present here the equivalent of figures 
131 and [5] to show that the results are essentially similar: 
as a reaction to the perturbation in the matter fluid, the 
DDE quickly forms a void. Starting very low, the rate of 



In this work we have investigated the clustering prop- 
erties of DDE. We modeled the DDE as a scalar field 
with a light mass, and have shown that in the vicinity of 
gravitationally collapsing matter, the DDE develops in- 
homogeneities and forms voids. Our results show a high 
sensitivity to the mass scale of the field. For a mass much 
larger or smaller than the Hubble scale, the field imitates 
the behavior of dust or the cosmological constant, respec- 
tively. The interesting dynamics is most prominent in the 
window where the mass and Hubble scales are compara- 
ble. This window is within the relevant mass range for 
dark energy models. As the Hubble function was larger 
in the past, heavier fields would have had comparable 
mass and Hubble scales at some point in the past. 

Our results should apply to any model of DDE which 
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achieves the present acceleration through a slow roll 
phase, as the slow roll assures that only a small patch 
of the potential is probed. We have shown this explicitly 
for the double exponential potential. Whether our results 
apply to an even wider class of models which achieve ac- 
celerated expansion not through slow roll should be fur- 
ther investigated. 

One thing which is clear from our results, is that DDE 
has potentially non-trivial behavior during the growth of 
inhomogencitics, though full non- linear analysis is needed 
to confirm whether the amplitude of the DDE inhomo- 
geneities is relevant to observations. As inhomogencitics 
of dark energy are a clear signature differentiating be- 
tween the cosmological constant and DDE, such possibil- 
ities should be fully explored and exhausted. 

A full treatment of the observational consequences of 
our results is beyond the scope of this work, but we would 
like to mention a few possibilities. Some obvious places to 
look for DDE inhomogeneities include lensing, the ISW 
effect in the CMB, number counts, and mass functions. 
Some of these directions are being pursued [1, Q . 

Our results show that both the energy density and the 
EOS of DDE develop a spatial dependence. Thus any 
observation constraining either of the above which can 
separately be measured locally and globally is valuable. 



Another possibility which we would like to point out, is 
that it would be useful to quantify the effect of a statisti- 
cal distribution of DDE voids on the CMB, or, following 
reference [l7|, on the directional distribution of super- 
novae. For example, detecting an angular inhomogeneity 
in H which is not in accordance with the matter distri- 
bution might suggest the presence of a DDE void. 
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